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NOTATIONS 


p 


Total  pressure  on  the  horizontal  strip  with  no  deflection 
Overpressure  on  the  soil 


Amount  of  arching 


Nondimensional  stress  = 


The  displacement  in  the  y  direction 

Critical  distance  from  the  center  of  the  strip  at  which 
the  net  pressure  on  the  strip  becomes  tensile 

Stresses  on  an  element  in  the  x  -  y  plane 


Strain  components  in  x  -y  plane 


Elastic  constant 


P 

P 

a ,  A. 

Y 

9 


Variables  of  integration 
Density  of  soil 
Airy's  stress  function 


INTRODUCTION 


If  the  foundation  of  an  underground  structure  settles,  the  pressure 
transmitted  to  the  structure  is  reduced  with  a  corresponding  increase  of 
pressure  in  the  neighboring  soil.  This  phenomenon  is  known  as  arching. 
Arching  due  to  deflection  of  a  rigid  horizontal  strip  or  base  has  been 
discussed  by  Terzaghi  ,  and  formulas  for  the  pressure  transmitted  to  the 
strip  are  given  based  on  assumed  failure  planes.  The  magnitude  of  the 
displacement  to  produce  these  failure  planes  is  not  known.  If  the 
deflection  of  the  strip  is  less  than  a  critical  value,  failure  planes 
will  not  develop,  and  the  arching  formulas  cannot  be  used  effectively. 

At  zero  displacement,  the  pressure  on  the  strip  is  equal  to  the  pressure 
of  the  soil  above  the  strip  plus  any  additional  overpressure  acting  on 
the  soil.  The  pressure  on  the  strip  decreases  as  the  displacement  is 
increased.  It  should  be  possible  to  displace  the  strip  to  a  critical 
value  such  that  all  the  pressure  acting  on  it  is  transferred  to  the 
neighboring  soil.  The  objectives  of  this  study  are  to  establish  the 
limits  of  the  critical  displacement  and  to  find  the  amount  of  arching  and 
the  configuration  of  the  pressure  distribution  on  the  base  as  a  function 
of  deflection. 
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W.  D.  Finn  has  treated  various  problems  dealing  with  stresses  in 
idealized  soil  media  subjected  to  different  types  of  boundary  conditions. 
One  such  problem  deals  with  the  stresses  in  soils  due  to  the  deflection 
of  a  rigid  horizontal  strip.  The  depth  of  the  soil  was  taken  as 
infinity,  which  imposes  a  restriction  in  adapting  the  solution  to  prac¬ 
tical  problems. 

The  present  study  deals  with  the  stresses  in  a  soil  field  of  finite 
depth,  h,  due  to  deflection  of  a  rigid  strip  of  width  2b.  The  soil  mass 
is  assumed  to  be  a  homogeneous,  elastic,  isotropic  medium  subjected  to 
high  overpressures.  If  there  is  no  deflection  of  the  horizontal  strip, 
the  pressure,  p,  transmitted  to  the  strip  will  be  equal  to  the  over¬ 
pressure,  p  ,  plus  the  pressure  of  the  soil  above  the  base,  Yh,  where  Y 
is  the  density  of  the  soil.  However,  if  the  strip  is  displaced  by  an 
amount  d;.  the  pressure  transmitted  to  it  will  be  reduced.  Assuming  that 
the  principle  of  superposition  is  valid,  the  amount  of  arching  (the 
amount  of  pressure  that  is  transferred  to  the  neighboring  soil  or  the 
reduction  of  pressure  on  the  strip)  is  equal  to  the  amount  of  tensile 
forces  on  the  base  due  to  the  displacement  d.  However,  as  the  strip 
deflects,  zones  of  very  high  tensile  stresses  form  toward  the  edges 
because  of  the  discontinuous  displacement.  Thus,  when  the  overpressure 
is  superimposed  on  the  tensile  stress  field,  there  will  still  be  residual 
censile  stresses  toward  the  edges  of  the  base.  Since  the  soil  media 
cannot  be  expected  to  transmit  tensile  stresses,  these  stresses  are  not 
considered  to  contribute  to  arching.  This  condition  is  specified  when 
computing  the  amount  of  arching. 


ANALYTICAL  DEVELOPMENT 


Procedure 


Figure  1  represents  a  section  through  a  soil  mass  of  depth  h  and 
infinite  width  subjected  to  an  overpressure  p  .  The  distance  in  the  z 
direction  perpendicular  to  the  x-y  plane  is  considered  as  infinity. 

The  width  of  the  rigid  horizontal  strip  is  taken  as  2b,  and  the  amount 
of  plate  deflection  given  by  d.  The  total  pressure  on  the  strip  with 
conditions  given  in  Case  (a),  Figure  1,  can  be  considered  to  be  equiva¬ 
lent  to  the  superposition  of  Cases  (b)  and  (c) .  Case  (b)  represents  a 
uniform  compression  on  the  base  exerted  by  the  overpressure  and  the  soil 
above  the  base  (p  =  p  +Yh),  with  no  deflection  of  the  rigid  strip. 

The  distribution  of  pressure  in  Case  (b)  is  indicated  by  line  12.  The 
tensile  stress  distribution  due  to  the  displacement  d  alone  is  assumed 
to  be  given  by  34567  of  Case  (c) .  The  tensile  stresses  at  the  edge  are 
infinite  due  to  the  discontinuity  of  displacement.  The  pressure  distri¬ 
bution  for  Case  (a)  can  be  assumed  to  be  given  by  the  superposition  of 
Cases  (b)  and  (c)  for  small  displacements  of  the  strip. 

It  is  to  be  noted  that  the  tensile  stresses  due  to  the  base 
displacement  reach  the  value  of  the  maximum  compression,  p,  at  some 
critical  distance,  +x  ,  from  the  center  of  the  base.  Beyond  this 
region  the  net  pressure  on  the  base  is  tensile.  These  resultant  tensile 
stresses  in  the  region  beyond  +x  from  the  center  of  the  base  are  not 
considered  effective;  thus,  the  net  compressive  force  acting  on  the 
strip  is  given  by  the  area  456. 

The  express  objectives  of  this  study  are  to  find  the  distribution 
of  tensile  stresses  as  shown  in  Case  (c)  by  34567,  to  determine  the 
distance  +x  at  which  the  resultant  pressures  become  tensile,  and  to 
find  the  amount  of  arching  as  shown  in  Case  (a)  by  14562. 

Evaluation  of  Stress  Function 


Since  it  is  assumed  that  there  are  no  strains  in  the  z  direction, 
the  problem  can  be  considered  as  one  of  "plane  strain,"  and  the 
appropriate  equations  of  the  theory  of  elasticity  can  be  used.~*  The 


positive  directions  of  the  stresses 
Figure  1.  Since  the  surface  is  free 
conditions  for  Case  (c)  at  y  =  h  are 

a  =0 

y 


r  =  0 

xy 


•  ,  cr  ,  and  T  are  3hown  in 

ol  applied  pressure,  the  boundary 
given  by 

(1) 

(2) 
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The  frictional  resistance  at  the  base,  y  =  0,  can  be  assumed  to 
prevent  any  elongation  in  the  x  direction  at  that  level.  Thus,  the 
strain  in  the  x  direction  at  y  =  0  is 


e 

x.,y  =  o 


Rcr 


=  0 


where 


S3  = 


and 


_  M*  (1+h) 
E 


(3) 


(4) 


E  is  the  modulus  of  elasticity  of  the  soil  media,  and  M-  is  Poisson's 
ratio . 


The  displacement,  v,  in  the  y  direction  at  the  boundary  y  =  0  is 
given  by 


v  =  -d  for  -b  £  x  <:  +b‘ 

(5) 

=  0  for  x<-b  and  x>+b 

The  displacement  at  the  boundary  can  be  expressed  in  integral  form 
by  using  the  Eourier  cosine  integral: 


v(x)  = 


TT 


cos  ax  da 


v(\)  cos  a\  d\ 


Thus,  the  displacement  at  y  =  0  as  defined  by  Equation  5  can  be 
obtained  as 


(6) 


v(x) 


to 

Jo 


sin  qb 
a 


cos  ax  da 


(7) 


where  a  and  \  are  variables  of  integration. 

The  displacement,  v,  can  be  expressed  in  terms  of  strain  e  in  the 
y  direction  as  ^ 


e 


y 


irv 

dy 


Ray  -  Pax 


(8) 


and  v  = 


e  dy  +  g(x) 


(9) 
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To  solve  for  the  various  stresses  in  the  media,  with  the  given 
boundary  conditions,  a  stress  function^  in  the  following  form  is 
assumed : 


r  i 

2 

a 

Jo  U 


Ae°^  4-  Bae0^  +  Ce  +  D^e  cos  Q?x  dx  (10) 


where  A,  B,  C,  and  D  are  constants. 

It  can  be  shown  that  the  above  stress  function,  cp  ,  satisfies  the 
biharmonic  equation 

4  d^cp  .... 

V  cp  -  f"  +  2  — TtV  +  _  4  ‘  °  (11) 

ox  ox  oy  Ay 

The  stresses  in  the  media  are  related  to  the  stress  function  by  the 
following  relations: 


O'  =  - 

X  dx2 


a  --O- 

y  dy2 


xy  dx  dy 

The  constants  A,  B,  C,  and  D  are  found  by  using  the  boundary  con¬ 
ditions  (1),  (2),  (3),  and  (7)  as  specified  above.  The  functions  g(x) 
in  Equation  9  can  be  shown  to  be  zero  after  appropriate  substitutions 
and  integrations  using  the  condition  that  v  =  0  for  h  — >  ® .  The 
relations  for  the  constants  A,  B,  C,  and  D  are  given  below: 


A  =  C  ( - 2o'h  -  1)  +  D  (-2a2h2)  e_2Q,h 


B  =  2C  +  D  (2ah  -  1)  e 


4 


-  2  p 


(15) 


Substitution  of  the  expressions  for  A,  B,  C,  and  D  in  Equation  10 
and  subsequent  simplification  will  result  in  an  expression  for  the  stress 
function.  The  stresses  anywhere  in  the  medium  can  be  found  by  using 
Equation  12.  However,  since  the  q  stresses  over  the  base  are  of  main 
interest,  only  these  are  evaluated^in  the  next  section. 

Pressure  Distribution  Across  the  Width  of  the  Strip 

The  expression  for  the  stress  cr  is  obtained  by  substituting  the 
stress  function  in  Equation  12  as  follows: 


5 


a 

y 


This  expression  for  O'  can  be  specialized  to  obtain  the  pressure 
distribution  across  the  horizontal  strip  by  integrating  and  then  sub¬ 
stituting  y  =  0  in  the  final  results.  However,  with  some  care,  y  *  0 
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can  be  substituted  in  some  terms  of  the  numerator  of  the  integrand 
before  performing  the  integration,  and  the  following  integral  is 


obtained 


for  cr 

y,o 


along  y  =  0: 


-  10P2  +  4j3p  -  2P2 

-  3P2  -  20P  +  P2 


C.dE  =  -  2  p  C_ 
4  i 


A  99 

4  tte  (-  30  -  20P  +  P  ) 

,  C^,  and  are  dimensionless  constants,  and  has  the  units  of 
pounds  per  inch. 

The  integration  of  the  expression  given  by  Equation  18  is  achieved 
by  expanding  the  denominator  in  a  series.  Thereafter,  all  but  the  first 
six  terms  of  this  alternating  series  are  neglected.  The  resulting 
expression  is  shown  below: 


7 


a  -  C.dE 
y,o  4 


,  -2ah  2,2  ,  _  x 

1  -  e  (40^0-  h  +  C2) 


+  e 


-4o/h 


-  e 


-6ah 


2  2  2 

(40^^  hz  +  c2r  -  i 


(4Clcv2  h2  +  C2)3  -  2  (4CjQ'2h2  +  c  ) 


4-  e 


-8ah 


2  2  4 

r  /-i  \  “ 


2.  2  .  „  v  2 


(40^  h  +  C2)  -  3  (40^  n  +  C2)  +  1 


-  e 


-10ah 


.2U2  ,  „  x5  _  /A„  „2,  2  .  „  s 3 


(40^^  +  C2y  -  4  (4C1«,th4  +  C2)J  +  3  (4Cia2h2  +  Q  ) 


+  e 


■12ah 


+  I  e  +  4ah  e 


“4o'h  .  -2*h  _  -ay, 


2  sin  ta  cos  da  (19) 


or 


a  =  S  a(l) 
y,°  y,o 


=  ct(1)  +a(2)  +a(3)  +a(4)  +a(5)  +  a<6>  + 
y,o  y  9  o  y  9o  y^o  y^o  y?o 


(20) 


where 


=  C,  dE 

r 

1 

y  > 0  4 

, 

0 

_  _ 

-4a'h  ,  -2o?h  -ay 

e  +  4ah  e  -  e  J 


2  sin  ba  cos  xa 


da 
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« 


(2)  _  ,F 
a  -  -C/L  uE 

y,o  4 

CO 

-2ah 

e 

2  2 

4C  aTi  +  C2 

-4oh  ^  ,  ,  -2ah  -oy 

e  +  4oh  e  -  e 

Vi 

0 

- 

—  — 

2  sin  bo  cos  xo 


do 


c/3)  =  C,  dE 
y,o  4 


-4oh 


2  2  2 
(4ClQ  h  +  C2) 


-4oh 

e 


+  4oh  e 


-2a'h 


2  sin  bo  cos  xo 


da 


- 

— 

a  (4)  =  -C,  dE 
y,o  4 

r 

-6oh 

e 

?  9  3 

+  c2)J 

-  2  (40^ 2h2  +  C2) 

0 

- 

- 

-4oh 

e 


+  4oh  e 


-20h 


-oy 


“ 

2  sin  bo  cos  xo 


do 


o(5)  = 

y,° 


C.  dE 
4 


-  8oh 
e 


•’ o 


2  2  4 

(40^  h  +  C2)4 


3  (4C1o2h2  +  C2)2  +  1 


-/Ph  ,  .  ,  -20h  -Oy 

e  +  4oh  e  -e 


[ 


2  sin  bo  cos  xo 


do 


9 


>co 


7<6>  = 
y,° 


-  C.  dE 
4 


-LOah 


?  2  5  2  2  3 

(4c  +  c2r  -  4  (40^ v  +  c2r 


+  3(4C  (A2  +  C2) 


[ 

[ 


-  4ah  .  -  2ah 

e  +  4ah  e 


-ay 


2  sin  bo?  cos  xa 


da 


,(i) 


(21) 


The  expressions  for  a  are  evaluated  individually  and  then  y  is 


y>° 


(1) 


set  to  zero.  In  the  expression  for  Q,  if  the  value  of  zero  for  y  is 
substituted  in  the  integrand,  the  improper  integral 
r*° 

sin  bo*  cos  xa  da 

4-  -ay 

does  not  exist.  However,  with  the  term  e  in  the  integrand,  the 
value  of  the  integral  is  given  by 


-  C.  dE 
4 


00 

■* 

e  (2  sin  be*  cos  xa)da  = 

J  o 


C.  dE 
4 


_ b  +  x _ 

2  2 
y  4-  (b  +  x) 


+ 


_ b  -  x _ 

2  2 
y  +  (b  -  x) 


and  the  limiting  value  as  y  0  is  given  by 


(22a) 


-  C.  dE 
4 


I 

b  +  x 


(22b) 
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Selected  tabulated  results  and  a  consideration  of  the  errors  of  the 
solutions  are  given  in  Appendix  A-II.  For  practical  purposes,  the 
accuracy  is  deemed  quite  adequate. 

The  results  showing  S  versus  x./b  for  different  values  of  b/h  and  p, 
are  given  in  Figures  2  through  5.  In  all  these  cases  the  value  of  S 
reaches  infinity  at  x/b  =  1.0  due  to  the  discontinuity  at  the  boundary. 
For  smaller  values  of  b/h,  the  curve  in  the  central  portion  if  flatter 
chan  for  larger  values  of  b/h.  for  the  extreme  case  of  b/h  :  1.0  and 
u  =  0.5,  the  pressure  at  the  center  of  the  base  (at  x/b  =  0)  is  very 
small  in  compression,  whexeas  for  other  portions  of  the  strip  the 
stresses  are  tensile.  Comparing  the  plots  for  various  values  of  u,  it 
can  be  seen  that  there  is  little  difference  in  the  stress  distribution 
except  for  the  extreme  cases  of  b/h  =  1.0. 
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Computation  of  Arching 

The  intensity  of  the  resultant  pressure  at  any  point  on  the  strip 
when  the  soil  of  depth  h  is  subjected  to  an  overpressure  of  p  is  given 
by  0 


p  , .  .  =  p  +  Yh  -  a  =  p  -  a 

^resultant  ro  y,o  y,o 

Nondimensionalizing  by  E(d/h),  Equation  25  can  be  rewritten  as 

(  £ resultant^  h  _  ph 

dE  dE  '  b 


(25) 


(26) 


The  term  Presuitant  Equations  25  and  26  should  be  considered  as 
equal  to  zero  for  negative  values  of  (ph/dE)  -  S. 


The  ratio  +x  /b,  for  which  the  resultant  pressure  becomes  negative, 
can  be  found  for  any  given  set  of  b/h,  p.,  and  ph/dE  parameters.  Thus, 
the  total  tensile  forces,  R,  over  the  half-width  of  strip  due  to  dis¬ 
placement,  d,  can  be  obtained  by  integrating  six  stress  terms  as  given 
by  Equation  30  (Appendix  A)  with  respect  to  x,  evaluating  the  definite 
integral  from  0  to  +  x^,  and  adding  the  tensile  forces  from  x^  to  b.. 

Thus,  the  amount  of  arching,  that  is  the  total  pressure  transferred  to 
the  neighboring  soil,  is  given  by 

R  =  2  P  (b  -  xcr)  (27) 


where 


Thus , 


percentage  of  arching  = 


R 

pb 


100 


(28) 


A  +  1 


x 

cr 

b 


100 


where 


A  = 


6 

£ 

i=l 


A. 

i 


6 

£ 

i=l 


pb 
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A  =  A  ft  !cr  £b 

i  i  h’  b  ’  ’  dE; 


The  expressions  for  AVs  are  given  in  Appendix  B-I. 


Given  the  parameters  ph/dE,  M-  ,}  and  b/h,  a  digital  computer  program 
was  developed  to  find  the  value  of  x  /b.  The  program  initially  assumes 

a  value  of  x/b  =  0,  computes  the  value  of  S,  and  compares  it  with  ph/dE. 

If  the  value  of  S  is  greater  than  ph/dE,  it  indicates  that  the  nqt 
pressure  on  the  base  is  tensile  and  an  arching  of  100  percent  is  indi¬ 
cated.  If  the  value  of  S  at  x/b  =  0  is  smaller  than  ph/dE,  then  a 
certain  increment  is  given  to  x/b  and  a  new  value  of  S  (x/b,  b/h, ^  )  is 
computed  and  is  compared  with  ph/dE.  If  the  difference  between  these  two 

is  less  than  or  equal  to  10  that  value  of  x/b  is  taken  as  x  /b.  If 

-5  cr 

the  difference  is  greater  than  10  ,  the  program  assumes  another  value  of 

x/b  and  the  process  is  repeated  until  the  value  of  x  /b  is  reached.  For 

cr 

values  of  x/b  approaching  1,  the  computer  takes  a  very  long  time  to  find 

the  value  of  x  /b.  In  cases  where  the  number  of  iterations  exceed  more 
cr 

than  100,  the  computer  prints  out  "number  of  iterations  more  than  100  to 
find  xcr"  and  proceeds  with  the  next  problem.  After  the  value  of  x^/b 

is  computed,  the  expressions  for  the  A^’s  are  evaluated  and  the  percentage 

of  arching  is  computed  using  Equation  28.  For  any  problem,  the  values  of 
b/h,  p  ,  ph/dE,  x^/b,  A-^,  A^,  A^ ,  A^,  A,-,  Ag,  A,  and  arching  are  printed. 

'The  following  values  are  used  for  the  parameter  ph/dE,  and  the  range  is 
considered  adequate. 


0.01 

0.1 

1 

10 

100 

0.0125 

0.125 

1,25 

12,5 

125 

0.015 

0.15 

1.5 

15 

150 

0.02 

0.2 

2 

20 

200 

0.03 

0.3 

3 

30 

300 

0.04 

0.4 

4 

40 

400 

0.05 

0,5 

5 

50 

500 

0.06 

0.6 

6 

60 

600 

0,07 

0.7 

7 

70 

700 

0.09 

0.9 

9 

90 

900 

1000 


Selected  tabulated  results  and  a  study  of  the  errors  involved  in 
taking  the  first  six  terms  of  the  series  are  given  in  Appendix  B-II. 
The  errors  for  the  most  part  are  found  to  be  small. 


The  value  of  x  /b  is  zero  when  the  arching  is  100  percent,  and  at 
zero  percent  arching  the  value  of  x  /t»  is  very  nearly  1.  Figures  6 
through  9  indicate  the  variation  of  x_^/b  with  ph/dE  for  Poisson's  ratio 
equal  to  0.1,  0.25,  0,3333,  and  0.5- 
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Figures  10  through  13  indicate  the  percentage  of  arching  versus  the 
yield  parameter  ph/dE  for  different  values  of  b/h  and  Poisson's  ratio, 

|i.  For  each  value  of  b/h,  the  percentage  of  arching  decreases  with 
increasing  values  of  ph/dE. 

Comparing  the  plots  of  arching  versus  ph/dE  and  x  /b  versus  ph/dE 

for  different  values  of  Poisson's  ratio,  it  can  be  seen  that  the  effect 
of  Poisson's  ratio  can  be  neglected  over  a  wide  range  of  ph/dE  for 
small  values  of  b/h. 

Illustrative  Example 

It  is  desired  to  find  the  amount  of  arching  developed  when  a  rigid 
horizontal  strip  24  feet  wide  buried  under  17  feet  of  soil  cover  under¬ 
goes  a  displacement  of  2  inches.  The  soil  is  subjected  to  an  overpressure 
of  100  psi.  The  modulus  of  elasticity  of  the  soil,  E,  is  10,000  psi, 
Poisson's  ratio,  p,,  is  0.25,  and  the  density  of  the  soil,Y  ,  is  110  pcf. 
Thus , 


b  =  12  ft 
h  =  17  ft 
b/h  =  0.706 
Pq  =  100  psi 


Y  =  no  pcf 


Yh 

p 


noa7j 

144 


13  psi 


p  +  Yh  =  113  psi 
o 


...  Ph  113  ('17Y  f  121  .  __ 

yreld  parameter,  ^  **  1>15 

From  the  plots  in  Figure  7  the  value  of  x  /b  =  0.77  and  from 

cr 

Figure  11  the  amount  of  arching  can  be  found  to  be  52  percent  for 

b/h  =  0.706  and  ph/dE  =  1.15.  The  distribution  of  pressure  on  the  base 

can  be  obtained  from  Figure  3.  At  x/b  *  0,  the  pressure,  p,  at  zero 

deflection  is  reduced  by  a  =20,6  psi,  where  u  is  obtained  from 

y,o  y,o 

S  =  crh/dE  =  0.21.  Thus,  the  net  pressure  on  the  base  at  x/b  =  0  is 
113  -  20.6  =  92.4  psi,  and  the  pressure  reduces  to  zero  at  a  distance 
+0,77  x  12  =  +  9.24  feet  from  the  center  of  the  base. 
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CONCLUSIONS 


The  analysis  indicates  that  arching  for  the  case  considered  here 
varies  from  100  percent  to  zero  percent,  depending  upon  the  three 
parameters,  b/h,  ph/dE,  and  (j, .  However,  for  practical  purposes  the 
effect  of  Poisson's  ratio,  ,  can  be  neglected  over  a  wide  range  of 
parameters.  The  errors  found  in  taking  the  first  six  terms  of  the 
infinite  series  are  small  for  a  majority  of  cases.  Experimental  data  to 
determine  the  validity  of  the  theory  are  not  available  at  this  time. 
However,  these  results  were  compared  with  the  data  obtained  for  a  case 
very  similar  to  the  one  considered  here,  and  the  trends  for  the  amount 
of  arching  are  similar  to  those  given  here. 
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Appendix  A 


Si  SERIES,  EQUATION  24 


I.  EXPRESSIONS  FOR  S. 's 

1 

The  following  notations  are  used  when  expressing  S_^'s: 


A(i)  (J)  = 


B(i)  (J)  = 


A(i)2(J)  = 


B(i)2(J)  = 


J’ 


J 


(29) 


The  last  two  notations  are  used  only  in  the  denominators  of  the  expres¬ 
sions  for  S..  's  and  i  kj-1.  The  resulting  relations  are: 


a  (1) 

s  =  =  c 

S1  F  d 

E  h 


0.25 


A41  B41_ 

A421  B421 


Al^  BU_ 
)  A222  B222 


1 


All  Bll 


(30a) 


16 


6(661)  -  20 (B63)  +  6(B65) 
B626 


-  0.375 


5 (A41)  -  10 (A43)  +  A45 
A425 


+  5(B41) 


10 (B43)  +  B45 
B425 


J 


1 


+  Clc2 


0.03125 


3  (A81)  -  A83 
A823 


+ 


3 (B81)  -  B83 
B823 


+  0.1481481 


4(A61)-  4(A63)  4(B61)  -  4(B63)  n 

'  A624  B624 


-  0.25 


3 (A41)  -  A43  3 (B41)  -  B43 

A423  B423 


] 


+ 


0.1111111 


2 (A61) 
A622 


+ 


2 (B61) 
B622 


-  0.25 


A41 

A421 


(30c) 


18 


(4) 


0.0046080 


7CA101)  -  35 CA103)  +  2KA105)  -  A107 
A1027 


7  ('B101')  -  35(B103)  +  21fBl05')  -  B107 
B1027 

J 


-  0.0769043 


8fA8:n  -  56(483")  +  56fA85^)  -  8(A87) 
A828 


8CB81-)  -  56 (B83')  +  56(B85)  -  8(B87) 
+  B828 


+  0.1646091 


7  (A6l~)  -  35  (A631  +  21(A65)  -  A67 
A627 


7(B61)  -  35fB631  +  2KB651  -  B67 
B627 


+  cl2c2 


0.01152 


5CA101)  -  10(A103)  +  A105 
A1025 


5fB10n  -  10 ('B103')  ±  B105 
B1025 


-  0.0878906 


6fA8n  -  20CA33)  +  6(A85) 
A826 


19 


20 


21 


22 


0) 

o 

CO 


* 


23 


nnnr1^c  r  1KA141)  -  165(Ai43)  +  462(A145)  -  330(A147)  +  55(A149)  -  A1411 

0.000917b  A14211 


1 1  ('B141')  -  165 fB143~)  +  462('B145')  -  330(B147)  +  55(B149)  -  B1411 


24 


26 


+  <c25  -  4c23  +  3<v  [_■  °-°714286  iitk +  iuii 


n  no77778  umn  ,  mm: 

0.0277778  A1222  +  B1222 


+  0  10  +  B101_ 

A1021  B1021 


(30f) 


II.  ACCURACY  OF  THE  SOLUTIONS 

A  digital  computer  evaluation  of  Equation  24  was  made  of  the  first 
six  S^  terms  for  the  parameters  b/h,  |i ,  and  x/b  to  study  the  possible 

errors  in  the  alternating  series.  Typical  results  are  given  in  Tables  I 
through  VII  for  several  values  of  b/h. 

It  is  extremely  difficult  to  prove  the  convergence  of  the  series. 

If  the  sixth  term  is  sufficiently  small  compared  to  the  first  term,  the 
sum  of  the  series  is  assumed  sufficiently  accurate  for  engineering  pur¬ 
poses  . 

For  b/h  =  1.0,  p,  =  0.1,  and  x/b  =  0,  the  value  of  is  0.1207  and 

Sg  is  -0.0310,  and  the  ratio  of  S^/S^  is  0.26.  A.lso  the  convergence 

appears  to  be  slow.  For  x/b  =  0.9,  the  value  of  S,  is  1.9630  and  S, 

i  o 

is  -0.0264,  and  the  ratio  of  js^/S^  is  0.013.  Even  though  the  sixth 

term,  S^,  is  considerably  smaller  than  the  first  term,  it  is  larger  in 

magnitude  than  S,_  by  a  small  amount.  However,  for  engineering  purposes 

the  value  of  S  =  1.9259  is  considered  adequate.  For  b/h  =  1.0,  p,  =  0.5, 
and  x/b  =  0,  the  ratio  of  |  S^/S^ |  is  0.35.  However,  observing  the  terms 

in  the  series,  it  can  be  seen  that  the  convergence  is  more  rapid  than  in 
the  case  where  p,  =  0.1.  As  the  value  of  b/h  decreases,  the  ratio  of 
Js^/S^  j  decreases  rapidly.  For  x/b  =  0,  the  values  of  |  S ^ / S ^ j  are  given 

for  various  parameters: 
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b/h 

jj.  =  0.1 

H  =  0.5 

0,05 

0.0002 

0.0003 

0.10 

0.0008 

0.0013 

0.20 

0.0034 

0.0053 

0.30 

0.0082 

0.0125 

0.40 

0.0155 

0.0237 

0.50 

0.0267 

0.0401 

0.60 

0.0430 

0.0641 

0.70 

0.0676 

0.0987 

0.80 

0.1043 

0.1417 

0.90 

0.1617 

0.2274 

1.00 

0.2568 

0.3521 

These  ratios  are  the  maximum  possible,  and  they  decrease  rapidly 
as  the  value  of  x/b  increases.  By  observing  the  values  of  S. 's  for 
the  four  tabulated  values  of  x/b,  it  can  be  seen  that  the  convergence 
is  slower  for  small  values  of  Poisson's  ratio,  |x .  Also  the  influence 
of  the  first  term  increases  rapidly  as  the  value  of  x/b  increases. 

Thus,  for  practical  purposes,  the  accuracy  obtained  by  taking  the  first 
six  terms  in  the  series  is  considered  quite  adequate. 
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Table  I.  Stress  Distribution  Across  the  Base  for  b/h  =.  0.05 


n 

x/b 

Si 

S2 

s3 

s4 

s5 

s6 

s 

0.1 

0 

7.9921 

-0.0143 

0.0042 

-0.0023 

0.0018 

-0.0017 

7.9798 

0.3 

8.7847 

-0.0143 

0.0042 

-0.0023 

0.0018 

-0.0017 

8.7724 

0.6 

12.4797 

-0.0143 

0.0042 

-0.0023 

0.0018 

-0.0017 

12.4874 

0.9 

42.1545 

-0.0142 

0.0042 

0.0023 

0.0018 

-0.0017 

42.1425 

0.25 

0 

7.6192 

-0.0152 

0.0042 

-0.0020 

0.0013 

-0.0009 

7.6065 

0.3 

8.3747 

-0.0152 

0.0042 

-0.0020 

0.0013 

-0.0009 

8.3621 

0.6 

11.9164 

-0.0151 

0.0042 

-0.0020 

0.0013 

-0.0009 

11.9038 

0.9 

40.1873 

-0.0151 

0.0042 

-0.0020 

0.0013 

-0.0009 

40.1748 

0.3333 

0 

7.6191 

-0.0168 

0.0048 

-0.0022 

0.0013 

-0.0009 

7.6053 

0.3 

8.3747 

-0.0168 

0.0048 

-0.0020 

0.0013 

-0.0009 

8.3609 

0.6 

11.9163 

-0.0168 

0.0048 

-0.0020 

0.0013 

-0.0009 

11.9025 

0.9 

40.1871 

-0.0168 

0.0048 

-0.0020 

0.0013 

-0.0009 

40.1734 

0.5 

0 

8.4657 

-0.0275 

0.0102 

-0.0057 

0.0038 

-0.0027 

8.4437 

0.3 

9.3052 

-0.0275 

0.0102 

-0.0057 

0.0038 

-o.oq27 

9.2833 

0.6 

13.2404 

-0.0275 

0.0102 

-0.0057 

0.0038 

-0.0027 

13.2185 

0.9 

44.6526 

-0.0274 

0.0102 

-0.0057 

0.0038 

-0.0027 

44.6307 
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Table  II.  Stress  Distribution  Across  the  Base  for  b/h  =0.1 


x/b 

S1 

s2 

S3 

S4 

S5 

S6 

S 

0.1 

0 

3.9643 

-0.0284 

0.0084 

-0.0046 

0.0035 

-0.0033 

3.9399 

0.3 

4.3606 

-0.0284 

0.0084 

-0.0046 

0.0035 

-0.0033 

4.3362 

0.6 

6.2183 

-0.0282 

0.0084 

-0.0046 

0.0035 

-0.0033 

6.1941 

0.9 

21.0460 

-0.0279 

0.0083 

-0.0046 

0.0035 

-0.0033 

21.0220 

0.25 

0 

3.7793 

-0.0302 

0.0085 

-0.0040 

0.0025 

-0.0019 

3.7542 

0.3 

4.1572 

-0.0301 

0 . 0084 

-0.0040 

0.0025 

-0.0019 

4.1321 

0.6 

5.9281 

-0.0299 

0.0084 

-0.0040 

0.0025 

-0.0019 

5.9033 

0.9 

20.0638 

-0.0295 

0.0084 

-0.0040 

0.0025 

-0.0019 

20.0393 

0.3333 

0 

3.7793 

-0.0335 

0.0096 

-0.0045 

0.0026 

-0.0017 

3.7518 

0.3 

4.1571 

-0.0334 

0.0096 

-0.0045 

0.0026 

-0.0017 

4.1297 

0.6 

5.9281 

-0.0332 

0.0096 

-0.0045 

0.0026 

-0.0017 

5.9009 

0.9 

20.0637 

-0.0328 

0.0095 

-0.0044 

0.0026 

-0.0017 

20.0369 

0.5 

0 

4.1993 

-0.0547 

0.0203 

-0.0113 

0.0075 

-0.0055 

4.1556 

0.3 

4.6191 

-0.0546 

0.0203 

-0.0113 

0.0075 

-0.0054 

4.5755 

0.6 

6.5868 

-0.0541 

0.0202 

-0.0113 

0.0075 

-0.0054 

6.5436 

0.9 

22.2932 

-0.0535 

0.0201 

-0.0113 

0.0075 

-0.0054 

22.2505 

30 


Table  III.  Stress  Distribution  Across  the  Base  for  b/h  =  0.2 


M- 

x/b 

S1 

s2 

S3 

S4 

S5 

s6 

S 

0.1 

0 

1.9198 

-0.0556 

0.0166 

-0.0092 

0.0070 

-0.0066 

1.8721 

0.3 

2.1183 

-0.0551 

0.0166 

-0.0092 

0.0070 

-0.0066 

2.0711 

0.6 

3 . 0483 

-0.0537 

0.0164 

-0.0091 

0.0070 

-0.0066 

3.0023 

0.9 

10.4641 

-0.0515 

0.0161 

-0.0090 

0.0069 

-0.0066 

10.4200 

0.25 

0 

1.8302 

-0.0588 

0.0167 

-0.0080 

0.0050 

-0.0037 

1.7814 

0.3 

2.0195 

-0.0583 

0.0167 

-0.0080 

0.0050 

-0.0037 

1.9712 

0.6 

2.9061 

-0.0567 

0.0165 

-0.0080 

0.0050 

-0.0037 

2.8591 

0.9 

9.9757 

-0.0542 

0.0162 

-0.0079 

0.0049 

-0.0037 

9.9311 

0.3333 

0 

1.8302 

-0.0652 

0.0190 

-0.0089 

0.0052 

-0.0035 

1.7768 

0.3 

2.0195 

-0.0646 

0.0189 

-0.0088 

0.0052 

-0.0035 

1.9666 

0.6 

2.9061 

-0.0628 

0.0187 

-0.0088 

0.0051 

-0.0034 

2.8548 

0.9 

9.9757 

-0.0599 

0.0183 

-0.0087 

0.0051 

-0.0034 

9.9270 

0.5 

0 

2.0336 

-0.1064 

0.0400 

-0.0225 

0.0149 

-0.0108 

1.9488 

0.3 

2.2439 

-0.1053 

0.0399 

-0.0224 

0.0149 

-0.0108 

2.1601 

0.6 

3.2290 

-0.1022 

0.0393 

-0.0222 

0.0148 

-0.0108 

3.1478 

0.9 

11.0841 

-0.0972 

0.0383 

-0.0218 

0.0146 

-0.0106 

11.0074 

Table  IV.  Stress  Distribution  Across  the  Base  for  b/h  =  0.4 


V- 

x./b 

S1 

s2 

s3 

s4 

S5 

S6 

S 

0.1 

0 

0.8436 

-0.1011 

0.0321 

-0.0180 

0.0139 

-0.0131 

0.7522 

0.3 

0.9455 

-0.0981 

0.0316 

-0.0179 

0.0138 

-0.0131 

0.8618 

0.6 

1.4182 

-0.0894 

0.0305 

-0.0175 

0.0136 

-0.0130 

1.3423 

0.9 

5.1380 

-0.0761 

0.0285 

-0.0169 

0.0133 

-0.0128 

5.0741 

0.25 

0 

0.8042 

-0.1063 

0.0320 

-0.0157 

0.0099 

-0.0073 

0.7168 

0.3 

0.9014 

-0.1029 

0.0316 

-0.0155 

0.0098 

-0.0073 

0.8170 

0.6 

1.3520 

-0.0932 

0.0302 

-0.0151 

0.0096 

-0.0072 

1.2762 

o 

VO 

4.8982 

-0,0783 

0.0280 

-0.0144 

0.0093 

-0.0071 

4.8358 

0.3333 

0 

0.8042 

-0.1175 

0.0362 

-0.0172 

0.0101 

-0.0068 

0.7091 

0.3 

0.9014 

-0.1136 

0.0356 

-0.0170 

0.0101 

-0.0068 

0.8096 

0.6 

1.3520 

-0.1025 

0.0339 

-0.0165 

0.0098 

-0.0067 

1.2701 

0.9 

4.8982 

-0.0854 

0.0312 

-0.0157 

0.0095 

-0.0065 

4.8313 

0.5 

0 

0.8936 

-0.1903 

0.0759 

-0.0433 

0.0290 

-0.0212 

0.7436 

0.3 

1.0015 

-0.1837 

0.0745 

-0.0428 

0.0287 

-0.0210 

0.8572 

0.6 

1.5022 

-0.1646 

0.0704 

-0.0412 

0.0279 

-0.0205 

1.3742 

0.9 

5.4425 

-0.1353 

0.0639 

-0.0385 

0.0265 

-0.0197 

5.3393 
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Table  V. 


Stress  Distribution  Across  the  Base  for  b/h  =0.6 


x/b 

S1 

s2 

s3 

s4 

s5 

s6 

S 

0.1 

0 

0.4507 

-0.1303 

0.0452 

-0.0261 

-0.0204 

-0.0194 

0.3405 

0.3 

0.5243 

-0.1230 

0.0440 

-0.0257 

0.0202 

-0.0193 

0.4205 

0.6 

0.8553 

-0.1026 

0.0405 

-0.0246 

0.0196 

-0.0189 

0.7693 

0.9 

3.3586 

-0.0730 

0.0351 

-0.0227 

0.0186 

-0.0183 

3.2983 

0.25 

0 

0.4297 

-0.1354 

0.0448 

-0.0225 

0.0144 

-0.0108 

0.3201 

0.3 

0.4998 

-0.1274 

0.0434 

-0.0221 

0.0142 

-0.0107 

0.3972 

0.6 

0.8154 

-0.1049 

0.0394 

-0.0208 

0.0137 

-0.0104 

0.7322 

0.9 

3.2019 

-0.0724 

0.0332 

-0.0188 

0.0128 

-0.0099 

3.1467 

0.3333 

0 

0.4297 

-0.1487 

0.0503 

-0.0246 

0.0147 

-0.0100 

0.3114 

0.3 

0.4998 

-0.1396 

0.0486 

-0.0241 

0.0145 

-0.0099 

0.3894 

0.6 

0.8154 

-0.1141 

0.0438 

-0.0225 

0.0138 

-0.0095 

0.7268 

0.9 

3.2018 

-0.0773 

0.0363 

-0.0200 

0.0127 

-0.0089 

3.1446 

0.5 

0 

0.4775 

-0.2381 

0.1041 

-0.0612 

0.0415 

-0.0306 

0.2932 

0.3 

0.5554 

-0.2226 

0.1001 

-0.0595 

0.0406 

-0.0300 

0.3839 

0.6 

0.9060 

-0.1794 

0.0886 

-0.0546 

0.0380 

-0.0288 

0.7700 

0.9 

3.5576 

-0.1173 

0.0709 

-0.0469 

0.0338 

-0.0258 

3.4723 

33 


Table  VI*  Stress  Distribution  Across  the  Base  for  b/h  =  0.8 


M- 

x/b 

S1 

S2 

S3 

Q 

“4 

S5 

S6 

S 

0.1 

0 

0.2434 

-0.1416 

0.0554 

-0.0333 

0.0264 

-0.0254 

0.1249 

0.3 

0.3057 

-0.1308 

0.0530 

-0.0324 

0.0259 

-0.0251 

0.1964 

0.6 

0.5740 

-0.1005 

0.0462 

-0.0299 

0.0246 

-0.0243 

0.4902 

0.9 

2.4800 

-0.0574 

0.0360 

-0.0260 

0.0226 

-0.0229 

2.4323 

0.25 

0 

0.2320 

-0.1448 

0.0541 

-0.0283 

0.0185 

-0.0140 

0.1174 

0.3 

0.2915 

-0.1332 

0.0514 

-0.0274 

0.0180 

-0.0138 

0.1865 

0.6 

0.5472 

-0.1006 

0.0438 

-0.0248 

0.0169 

-0.0131 

0.4694 

0.9 

2.3643 

-0.0543 

0.0326 

-0.0208 

0.0150 

-0.0121 

2.3246 

0.3333 

0 

0.2320 

-0.1575 

0.0602 

-0.0307 

0.0187 

-0.0128 

0.1099 

0„3 

0.2915 

-0.1445 

0.0570 

-0.0296 

0.0182 

-0.0126 

0.1800 

0.6 

0.5472 

-0. 1080 

0.0478 

-0.0263 

0.0167 

-0.0118 

0.4657 

0.9 

2.3643 

-0.0563 

0.0344 

-0.0214 

0.0145 

-0.0106 

2.3249 

0.5 

0 

0.2578 

-0.2478 

0.1224 

-0.0749 

0.0518 

-0.0386 

0.0707 

0.3 

0.3239 

-0.2262 

0.1148 

-0.0714 

0.0499 

-0.0374 

0.1535 

0.6 

0.6081 

-0.1657 

0.0936 

-0.0616 

0.0444 

-0.0340 

0.4848 

0.9 

2.6270 

-0.0801 

0.0627 

-0.0468 

0.0359 

-0.0286 

2.5702 
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Table  VII.  Stress  Distribution  Across  the  Base  for  b/h  =1.0 


V- 

x/b 

S1 

S2 

S3 

S4 

S5 

S6 

S 

0.1 

0 

0.1207 

-0.1378 

0.0623 

-0.0392 

0.0318 

-0.0310 

0.0068 

0.3 

0.1771 

-0.1263 

0.0586 

-0.0377 

0.0310 

-0.0305 

0.0721 

0.6 

0.4106 

-0.0924 

0.0481 

-0.0335 

0.0287 

-0.0289 

0.3325 

0.9 

1.9630 

-0.0418 

0.0330 

-0.0271 

0.0251 

-0.0264 

1.9259 

0.25 

0 

0.1150 

-0.1380 

0.0597 

-0.0329 

0.0220 

-0.0159 

0.0089 

0.3 

0.1688 

-0.1261 

0.0556 

-0.0314 

0.0212 

-0.0165 

0.0717 

0.6 

0.3914 

-0.0909 

0.0443 

-0.0270 

0.0191 

-0.0153 

0.3217 

0.9 

1.8714 

-0.0377 

0.0282 

-0.0205 

0.0159 

-0.0135 

1.8438 

0.3333 

0 

0.1150 

-0.1481 

0.0656 

-0.0352 

0.0220 

-0.0153 

0.0040 

0.3 

0.1688 

-0.1352 

0.0608 

-0.0333 

0.0211 

-0.0149 

0.0673 

0.6 

0.3914 

-0.0965 

0.0474 

-0.0281 

0.0186 

-0.0135 

0.3194 

0.9 

1.8714 

-0.0378 

0.0285 

-0.0203 

-0.0148 

0.0114 

1.8452 

0.5 

0 

0.1278 

-0.2274 

0.1302 

-0.0838 

0.0594 

-0.0450 

-0.0387 

0.3 

0.1876 

-0.2068 

0.1193 

-0.0783 

0.0562 

-0.0429 

0.0351 

0.6 

0.4349 

-0.1448 

0.0892 

-0.0628 

0.0471 

-0.0370 

0.3266 

0.9 

2.0794 

-0.0500 

0.0470 

-0.0403 

0.0335 

-0.0280 

2.0416 

35 


Appendix  B 


A  SERIES,  EQUATION  27 


I.  EXPRESSIONS  FOR  A. 's 

1 

The  integrations  of  Equation  27  become  less  involved  if  they  are 
first  performed  with  respect  to  x  and  then  with  respect  to  O',  starting 
from  Equation  21.  The  expressions  for  A^'s  are  given  by  Equation  32. 

The  following  notations  are  used  in  Equation  32: 


C(i)(J) 


D(i)(J) 


C(i)2(J)  = 


where  i  i  j  and  the  last  two  notations  are  used  only  in  the  denominators 
except  when  used  with  logarithms. 


36 


37 


38 


39 


36(0102)  -  126  (C104)  +  84CC106)  -  9(C108 


J' 


0.0138889 


I 

I 


41 


D14210 


42 


0.0400091 


1  +  15  ('D102')  -  15  CD104)  +  D106 


II.  ACCURACY  OF  THE  SOLUTIONS 


Typical  results  are  given  in  Tables  VIII  through  XIV  for  the  digital 
computer  evaluation  of  Equation  28  for  the  first  six  A.  terms,  consider¬ 
ing  the  parameters  b/h,  pt 3  ph/dE,  and  x^/b.  As  for  tfte  S^'s,  the 

convergence  of  the  series  is  very  difficult  to  prove.  It  is  presumed 
that  the  results  are  acceptable  if  the  last  term  of  the  series  is 
sufficiently  small  compared  to  the  first  term.  From  the  results  it  can 
be  observed  that  the  error  is  large  for  large  values  of  b/h,  small 
values  of  ph/dE,  and  large  values  of  p,.  However,  the  convergence  is 
poor  for  small  values  of  pi.  The  following  table  represents  the  maximum 
ratio  of  jA^/A^j  for  all  values  of  b/h  where  q  =  0.5. 


b/h 

ph/dE 

Values  of  |A( 

0.05 

9.0 

0.0004 

0.1 

5.0 

0.0011 

0.2 

2.0 

0.0051 

0.3 

1.25 

0.0121 

0.4 

0.9 

0.0223 

0.5 

0.5 

0.0395 

0.6 

0.3 

0.0635 

0.7 

0.2 

0.0956 

0.8 

0.06 

0.1465 

0.9 

0.01 

0.2260 

1.0 

0.01 

0.3166 

For  other  values  of  ph/dE  and  pt  the  error  is  smaller  than  that  indicated 
above . 


Table  \zIIT.  Arching  Due  to  Deflection  of  Strip  for  b/h  =  0.05 
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Table  X..  Arching  Due  to  Deflection  of  Strip 
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Figure  6.  Distance  to  point  of  zero  pressure  on  the  base  for  fx  -  0.1. 


Figure  8.  Distance  to  point  of  zero  pressure  on  the  base  for  /i  =  0.3333. 


Figure 


« 


Values  of  ph/dE 

Figure  10.  Arching  due  to  deflection  of  the  base  for  (i  =  0.1. 


Values  olf  ph/dE 

Figure  11.  Arching  due  to  deflection  of  the  base  for  n  =  0.25. 


Values  of  ph/dE 

Figure  12.  Arching  due  to  deflection  of  the  base  for  fi  =  0.3333. 


Figure  13.  Arching  due  to  deflection  e>f  the  base  for  -  0.5. 
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